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ABSTRACT 


A computer-controlled data collection system was designed and 
implemented to analyze finite amplitude acoustic standing waves in a 
rectangular air-filled cavity with linear and wedge-shaped boundary 
perturbations. The response of a nonlinear standing wave in per- 
turbed and unperturbed cavities was studied. The experimental 
results showed that the computer-obtained data was in qualitative 
agreement with that previously obtained by manual procedures. Quan- 
titative comparisons with theory could not be made because the prop- 
erties of the cavity did not match those for which theoretical 


predictions are available. 
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I. INTR TION 


A one-dimensional mathematical model for finite-amplitude 
standing waves in a real cavity was studied by Coppens and Sanders 
[Refs. 1, 2] in 1968. In 1974, this model was extended to three 
dimensions. Losses and resonance frequencies, inadequately 
described by the classical theory of wall losses, are obtained for the 
real cavity from measurements obtained in the linear acoustic regime. 
The research done in 1974 showed that when a cavity is driven near 
the resonance of a standing wave the nonlinear distortion primarily 
involves only the normal modes of the cavity that are “familially 
related” to the driven mode. Lane [Ref. 3], DeVall [Ref. 4], Slocum 
[Ref.5], and Kuntsal [Ref. 6] examined the problem and showed that 
experimental results are in excellent agreement with the theory 
except when a degeneracy exists in the cavity. 

The purpose of the research reported in this thesis was to design 
and implement a computer-assisted data collection and analysis sys- 
tem for measuring the finite-amplitude distortion in a perturbed 


cavity with degeneracies. 
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II. BACKGROUND AND THEORY 


When the sound intensity in an acoustic process becomes so high 
that the principle of superposition no longer applies, we enter the 
domain of nonlinear acoustics. The limits of linear approximations in 
the equation of force, mass continuity, and state are exceeded and the 
resultant vibrational disturbance at any point in the medium will no 
longer be equal to the sum of its individual components. The essence 
of the nonlinear acoustic mechanism is that an intense wave perturbs 
the medium in which it exists and this perturbation alters the natural 
order governing the wave’s behavior. The result is that segments of 
the disturbance in a warmer compression travel faster than those in a 
cooler rarefaction. When this happens, the acoustical wave alters its 
Shape, steepening as it travels along. The alteration is infinitesimal for 
weak disturbances and gradually increases in importance with 
increases in both amplitude and frequency of sound. [Ref. 2] 

In 1968, Kirchoff [Ref. 7] studied the distortion of intense acous- 
tic waves in an unbounded fluid. This study was extended by Lamb 
(Ref. 8], Fay [Ref. 9], and by Keller [Ref. 10]. 

Keck and Beyer [Ref. 11] approached the problem of unbounded 
propagation by using perturbation methods and Fourier series 
representations of the waveform. The model of Coppens and Sanders 
[Ref. 12] dealt with finite-amplitude standing waves in rigid walled 
cavities. They extended the Keck-Beyer perturbation approach to 
include wall losses predicted by the Rayleigh-Kirchoff theory [Ref. 7], 


iat 


and showed excellent agreement with experiments conducted in a 
rigid-walled tube at low but finite levels of excitation. The 
experimental work by Beech [Ref. 13] and by Ruff [Ref. 14] showed 
that at high excitation levels differences developed between theory 
and experiment. Experimental work by DeVall [Ref. 4] and by Kilmer 
[Ref. 15] showed that if degeneracies existed the model failed to 
account for the experimentally observed excitation of “nonfamilial” 
modes. 

The effects of various geometrical boundary perturbations on a 
rigid-walled cavity investigated by Aydin [Ref. 16] showed that the 
standing waves that exist in an ideal cavity must be corrected when 
the boundaries are irregular. The model of Coppens and Sanders is 
based on a three-dimensional, nonlinear wave equation with a 
dissipative term describing absorptive losses encountered by standing 
waves in rigid walled cavities. As stated in [Ref. 12], 

The standing waves which can be excited in any real cavity deviate 
from predictions of the linear wave equation with ideal boundary 


conditions for the following reasons: 


a The presence of boundary-layer losses at the cavity surfaces 
yields a dispersive contribution to the wave equation. 


b. Geometrical irregularities alter the effective dimensions of the 
cavity. 


Both of these mechanisms can be treated separately as long as the 
shift in frequency is so small that the actual resonances are close to 
the theoretical values resulting from the classical model. These are 
treated by assuming the dimensions are exact, and the apparent phase 


speeds are determined on that basis. 


ie2 


The resonance frequency for each standing wave is defined as 


(Ref. 6] 


On = Cn [(nx kx)2 + (ny ky)? + (nz kz)?2]’? (1) 


l 
where Kx = _ ky = T kz, = = and n, m, and 1 are integers. Cn is 


the apparent phase speed appropriate for that frequency. 
The force equation, in Eulerian coordinates, can be written in the 


form 


a -V)a + AW (2) 


where the operator A describes those physical processes leading to 
absorption and dispersion. If we ignore rotational effects so that 
u = Vo, restrict waves to rectangular cavities, and neglect terms of 
orders higher than M2, M(a/k), and (a/k)2, the result is a 


quadratically nonlinear wave equation [Ref. 12]. 
ae P 02 |y-1f P r u 
Cp gai | rot 2 Ce. 3 (& | >) 


The standing pressure wave in a rectangular cavity has the form 


(Ref. 12] 


p = coskx x coskyy coskzz cost (4) 
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If the cavity is being driven near resonance, the nonresonant terms 
contribute almost nothing with respect to the resonant terms. So the 


standing wave pressure has the form of 
p= 2d Pn (5) 


where each Py is of frequency nw and has the form [Ref. 12]. 


Deed * MRy cosnkxx cosnkyy cosnkzz sin(nwt + Qn) (6) 


Po 


These pns constitute those normal modes familially related to p}. 


The nonlinear, coupled, transcendental equations which must be 


solved for the Ryps are 


om bid (Qn - 8n) = NMBQn cosa | by Rj Rag | SOS * * (0) — Gay 


=) Rnaj Ry {°° ost Ong 0) for alln > 1 (7) 


R=1l for n=] 


where Qr = Arctan(-F,) 


Fn=Q [(nw)2 - (on)? ]/( (nw)? = 2Qn a) for Ss << 1 (8) 
On 
Qn= Soop (9) 


Total pressure near resonance for perturbed boundaries is [Ref. 6] 
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= (n, m, 1/,0) ->—=—_ © am An (-1)" * 
on 
O 


Oo O 


© 


Snm! (n,m,1/@,8 + Onm)) (10) 
where Snm1 ~ first-order perturbation correction and Onm angle 


which is defined as [Ref. 6]: 


cs 


“dia, A = Amp. of perturbation 


If the cavity has perturbed walls, the solution will be in terms of a 
summation of the classical linear solution for ideal boundaries plus 


perturbation correction terms due to the irregular boundary: 


P=Pote p't+e2p"+... (11) 


Since the magnitude of the boundary perturbation is kept small, sec- 


ond and higher terms in e€ can be considered insignificant, so that 


p=Po+te p (12) 
and p must Satisfy the following conditions, 


Sey 0 (13) 
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and 


Wisi ot — (0, at x=0, Lx [1 + efly,z)] (14) 
y = 0, ly 
Z= Oa, 


where Ti, the local normal to the real surface, is obtained by taking the 
gradient of the equation for the boundary [Ref. 1], 


T=V &k-Lx [1 + efly.z)]} (15) 


Thus, to the first order in e, 


a8 et, MUAy_ op, Myals 16) 


and when equation (16) is used in equation (15), the result is 


ane 2te! oo a 


dz el Lileaeival 


Applying a Taylor-series expansion for p evaluated at the real 


boundary Lx [1 + efly.z)] produces 


= oP 
Pl Lx {1 + efly.z)] “x P | =~ * OX | 1, © lly.) 
,i% | le Lxfly.z)]2 + ... (18) 
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Substituting equation (12) into the RHS of equation (18), taking 
the partial derivative with respect to x on both sides and keeping the 


first-order terms in e, yields [Ref. 1] 


dp _ Po op' 02Do 
oma Ly {1 + efly.z)] = sel Gs +e Olt Oe ay, eat) + ... (19) 


Taking the partial derivatives with respect to y and z and using 


exactly the same procedure gives [Ref. 1] 








dp OPo (25 Do 

~ Ov . f(y, ee 
_ eee ee1.* a Le * OyOX, = | ees) (20) 
a = Pe | +e 2 +S Po | L, lata) + ... (ayy 
02x = Ly [1 + €fly.z)] OZ x=Lx az. =Ly * 9Z0X 


Substituting equations (19), (20), and (21) into equation (17) and 


keeping the first-order terms in e€ results in 


dp' - 22Po | afly.t) Po | afly.z) apa] 
Kye tg tly.2) Sa dy oy oz oz \x-L (22) 


The RHS of equation (22) can be represented as a Fourier series in 


cosines, so that p' can be expressed as a Summation of normal modes. 


Melee, 
ay = 3 * aml COS = y cos <= z cos (wt + 8) (23) 
sem on Eo Ly Ly 


Or 


1, 


a = py Ss am] (0,m,1/w,8) (24) 


where 
on el Sid dy dz (25) 
omo= 2) f 1c cos Ty dy Gh (26) 
2or= Bf f tcl cos Bay a a7 
aml Ac!" fia cos Ty cos i? dvds, Bae (28) 


and G is defined as 


£9 
sae ee (29) 


In order to find the contribution to p' from each of the terms, it is 
stated [Ref. 16] that, for a cavity forced by a dynamic boundary condi- 
tion at a boundary, 


Cae 
C UL p =0 (30) 


and the dynamic boundary condition from the (m,])th term is 
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dp’ IE 
al Acs ty cos 5 eilot+Q) (S35) 
then 
. a eae An (-1)2® S cos +x cos =~ peste 
P ml = a a, n nml ie oy, y ic. 
Co 
. ei(ot+Q+onm)) (32) 
where 


Applying this solution to equation (24) gives the complete solution 
for the first-order perturbation [Ref. 6] 





i [48 


aml An (-1)"® Snm 


De= "(i ae) — 
= L,. 10 
og m=o 
10 
(n,m,1/w,80+0nm}) (S383) 


As it is stated [Ref. 16], if this is near a resonance, ®~@pm}, then 
this term will dominate the summation and all other non-degenerate 


terms can be omitted. Consequently, equation (33) becomes 


p = (n,m,1/q@,8) en An (-1)® Spm) (n,m,1/@,90+0nm)) ... (34) 
Co Lx 
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Ill. EXPERIMENTAL NSIDERATION 


A. COMPUTER-CONTROLLED APPARATUS 

The rectangular cavity used in this research was designed and 
constructed in 1977 by Kuntsal [Ref. 6] (Figure 1). One wall of the 
cavity could be moved and secured. The interior of the cavity was 
12.002 inches long and 2.502 inches high, while the width could be 
varied from approximately 5.50 inches to 7.00 inches, in 0.25-inch 
steps. The piston, located in a 2.25-inch diameter port in the floor of 
the cavity, was supported by an O-ring to produce a seal between the 
piston and the port. Three 1/4-inch diameter microphone ports 
were used. A condenser microphone was mounted in a 0.24-inch 
aluminum case and an O-ring and silicon grease were used between 
this case and the ports for sealing. 

Modal designation for the lowest three modes is shown in 
Figure 2. The (020) and (100) modes are nearly degenerate. The 
(020) mode is a one-dimensional standing wave with nodes 3 and 9 
inches along the long wall of the cavity and the (100) mode is a 
standing wave with one node in the middle of the short wall. When 
the microphone is placed at position A, the (100) mode is sensed 
without contamination from the (020) mode. In position C, the (020) 
mode is sensed without contamination from the (100) mode. A block 
diagram of the apparatus is shown in Figure 3. 

The piston was driven by an M-B Electronics model EA1500 


exciter, which in turn was driven by an APS Electronic model 114 
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Microphone Port A 


Microphone 
Sieh Port B 
a H/ Microphone 








Port C 


Figure 1 


Rectangular Cavity 
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piston 
position 


(0,1,0) 


Mic. 


HP 87 
Computer 


B 
€ 
12" 
(a) 
(0,2.0) (1,0,0) 
(b) 
Prolines2 


(a) Cavity Orientation 
(b) Modal Designation 


Mic. 


simp. cos 5060A ree 
Multi. 
3478A 


Dyn.-Sig. Anl. 
HP 3562A 

I 

| 

3 ref. 

Function 

| Generator 
| 3325 Synt. 
| 


Figure 3 
Block Diagram 
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power amplifier. The driving signal was produced by an HP 3325A 
Synthesizer/Function generator. 

The piston movement was continuously analyzed by an Endevco 
model 2215 Accelerometer mounted within the piston. The output 
voltage of the microphone and output voltage of the accelerometer 
were measured by an HP 3678 multimeter. An HP 3562A Dynamic 
Signal analyzer was used to measure the harmonic distortion. A model 
59204 Lock-in Analyzer was used to measure the resonance frequen- 
cies at small signal amplitude. 

The sound pressure in the cavity was sensed by a Bruel and Kjaer 
4136 condenser microphone with matching preamplifier 2807. The 
output of this preamplifier went to three pieces of equipment: (1) an 
HP 3478 multimeter to measure the voltage level, (2) an HP 3562A 
Dynamic Signal Analyzer to measure the harmonics, and (3) a Cos 
9060A oscilloscope to monitor the undesired distortion of the accel- 
erometer signal. 

The HP 3562A Dynamic Signal Analyzer and HP 3325A Synthe- 
sizer/Function generator were controlled by an HP 87 computer. Two 
types of geometrical perturbations were used: 


a The long wall was machined at a small angle (as shown in 
Figure 4). 


b. The long wall was machined as shown in Figure 5. 
The effect of these geometrical perturbations on the resonance 
frequency of the (100) mode could be estimated by calculating an 


“effective” width Le = L - AL of the perturbed cavity where AL = V/A 
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Figure 4 
Linearly Perturbed Wall 





Figure 5 


Wedged Perturbation 


and V and A are the volume and cross-sectional area of the 


perturbation. [Ref. 1] 


B. STRENGTH PARAMETER (STRM) 


Strength of the finite amplitude interaction is characterized by 


the strength parameter STRM [Ref. 6]. 
STRM = MBQ, where Q is the quality factor of the driven mode 


| Cp 
= [1+ @)- 1.2 for air 


= oCz 
v, V2 


B 
M 





where Sm is the microphone sensitivity. 
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The microphone output voltage v,; is read on an HP 3562A 
Dynamic Signal Analyzer set at the fundamental harmonic. Micro- 
phone sensitivity was obtained with a Bruel and Kjaer model 4220 


piston phone. 


Sy leoo = 0.03) 1059 volt/(N/M-) 


By using p = 1.293 kg/m? and C = 345 m/sec, we found 


Sinn 7.079. 10-2 ay x OF [o)) 
This is the same sensitivity as is measured by Kuntsal [Ref. 6]. 


C. FREQUENCY PARAMETER (FP) 


The frequency parameter is defined by 


BP = 26) (@8] 1) / 1 (36) 


where f is a driving resonance frequency and f; is a resonance fre- 


Wn 





quency of (010) mode. It is assumed that aloe P imedicates 
where the driving frequency is with respect to the fundamental reso- 
nance frequency of the cavity. As |FP| increases, the nonlinear effects 


will in general become smaller for a given strength parameter. 
D. HARMONICITY COEFFICIENT (E(n)) 
The harmonicity coefficient is defined as [Ref. 6] 


n —~Nn@] 


as (37) 


DS) 


E(n) characterizes how well the modes of a given family are tuned. 


E. COMPUTER PROGRAMS 

Two computer programs were developed. The first program finds 
the resonance frequencies and calculates the quality factors (Figure 6). 
The second program drives the Dynamic Signal Analyzer HP 3562A by 
a HP 3325A Synthesizer/Function generator and finds and calculates 
harmonic content of the microphone output at each driving frequency 


(Figure 7). 
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Set freq « FR 










Read “FRI” 
Read “FR2” 
Read “FR3" 


Ke Kant 4 o> <4 00 


GOSUB BINDM 


FR <— MAX(FR) 


OSUB HALFPOWER 


L_ | <—— AMP HALFPOWER (FR) 


[AMP (MAX(FR)) + = 





CALCULATE °Q° 


Print MAX(FR) 
Print Halfpower (FR) 
Print (Q) 


Figure 6 
“Q” Calculation Flow Diagram 


Pre-set HP 3235A 


False 
560<FR<570 


N@ — POWER (FREQ) 


READ Fl < FREQ 







1. Harmonic Power - Nl 
2. Harmonic Power 7—- N2 
3. Harmonic Power + N3 









Read Nl 
Read N2 
Read N3 


E> 


GOSUB BINOM 
Find Nl — N@G 


Fl Gal +1025 


Read Nl Print NI, Fl 

Read N2 Print N2, F2 

Read N3 Print N3, F3 
Migure 7 


Power of the Harmonic Content Calculation Flow Diagram 
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IV. DATA LLECTION PROCEDURE 


A. SYSTEM ALIGNMENT 

The primary factor in setting the lower limit to the measurement 
of harmonic content was the amount of distortion injected into the 
cavity by the piston. Misalignment of the piston within its port causes 
distortion of the piston motion which was detected by analyzing the 
harmonic content of the accelerometer output voltage. (With the pis- 
ton removed from its port, distortion levels were of the same order of 
magnitude as the general electronic noise level of the accelerometer 
system when at rest.) Before each run, the piston was aligned within 
its port so that the accelerometer output voltage contained less than 
O.1 percent second harmonic and 0.2 percent third harmonic. 

After testing various configurations, the one selected for ease of 
alignment consisted of the cavity mounted on a tripod support directly 
above the exciter, which was supported by its own base in an upright 
position. Both the cavity tripod and the exciter rested on a 5/16-inch 
pad of stiff acoustic rubber. Once the system was aligned, only an 
occasional minor adjustment was necessary to maintain low piston 


distortion. 


B. AUTOMATIC DATA COLLECTION 

Since the resonance frequencies of the cavity were observed to 
drift with time, the system was allowed to warm up for at least three 
hours prior to data collection. The harmonic content of the 


accelerometer output was analyzed and the piston was adjusted until 
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the second harmonic of the accelerometer output was at least 50 dB 
below that of the fundamental frequency. Figure 8 shows the drift in f, 
is approximately linear in time. 

Data was collected using two measurement procedures in the 
computer programs: (1) a pre- and post-run infinitesimal amplitude 
measurement; and (2) the finite amplitude measurement. 


1. Pre-Run and Post-Run Infinitesimal Amplitude 
Measurements 


During the pre-run and post-run infinitesimal amplitude 
measurement, the piston was driven 80 mV (rms). Because the (100) 
and (020) modes were nearly degenerate, measurements of the prop- 
erties of the (100) mode were made with the microphone in port C 
while measurements of the (020) mode were made with the micro- 
phone in port B. The resonant frequency of the fundamental and the 


Q factor first few (0,n,0O) in the modes were determined from 


fu -f 
f= “3— (38) 


fr 
Q at ae (39) 


and the harmonicity coefficients were then calculated from 


fn - nf) 
gn — nh (40) 
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Frequency (Hz) 


969.2 eae?) 
ee pt mE 
= 


964,2 
064, | 





Time (ks) 


Figure 8 


Drift in Frequency 
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565 5.072 


10 


where f,, and f, are the frequencies of the upper and lower half-power 
points, and fp is the resonance frequency of the nth member of the 
(010) family. Figure 9 shows the precision of frequency 
measurements. 
2. The Finite-Amplitude Measurement 

All finite-amplitude measurements were made at port B. 
Throughout this measurement, V, had to be kept constant by the 
computer program (Appendix C), which required adjustment of the 
driving voltage applied to the piston. The frequency was increased by 
the computer in 0.5-Hz steps through +5 Hz about the resonance 
frequency of a fundamental mode. At each driving frequency, the 
harmonic content of the microphone output, usually up to the fourth, 
was measured with the HP 3562A Dynamic Signal analyzer. 

To determine the exact value of the frequency parameter, the 
value of fp at the instant of the measurement of Vyz was determined 


from the first computer program. 
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1 


964.46 


964.42 


964.38 





Figure 9 


Precision of Frequency Measurements 
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V. RESULTS 


In this section, the theoretical predictions and experimental 
results obtained by Kuntsal [Ref. 6] will be qualitatively compared to 
the experimental results obtained by using the computer-controlled 
equipment. 

The information on the empirical losses and resonance frequen- 
cies is contained in the Q’s and E’s. It is important to note that, 
because V is always kept constant as a function of frequency, the sec- 
ond harmonic distortion is maximized for o = @2/2. The point where 
this occurs for each Vo/V, curve is indicated by an arrow with the 


label Fo99. At that point, the value of the frequency parameter is 


Fo20 = 2 Qo10 Eo20 


Similarly, 


Fo30 = 2 Qo10 Eo30 


is indicated on each graph. The arrow labelled Fi0909 indicates the 


position where the nearly degenerate (1,0,0) mode is resonant, and 


the value of Figo is 


Fio0 = 2 Qo10 E100 


The theoretical predictions for unperturbed cavities made by 


Kuntsal [Ref. 6] were compared to the experimental results as seen in 
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Figure A. It is clearly seen that the theory and experiment are qualita- 
tively in agreement. 

At Figure 9, fj and fg are plotted as a function of time. This figure 
showed two things: (1) changes of resonance frequencies of the family 
modes are proportional to the frequency, as would be expected if the 
long dimension of the cavity changes with time, most likely due to 
temperature changes: and (2) the fluctuations of resonance frequen- 
cies from measurement to measurement indicate a precision of + 0.02 
Hz in determining the resonance frequencies of the lowest mode. By 
Memeo tie value a= ©-O2Z)hi2Zeexpected uncertainties in the E’s are 


f, + 0.02 ~- n(f; + 0.02) 


AE expected es nf 


Calculated AEexpecteq and measured AE's between pre- and post-run 
are tabulated in Table 4. The difference of the pre-run and post-run 
E’s for the family modes in the unperturbed cavity are within the 
expected precision. For the 100 mode (non-family mode), this differ- 
ence is much greater than expected, indicating that the cavity may be 
changing its dimensions unequally in all directions. 

For a linear perturbation (Figure B) and for the wedge perturba- 
tion (Figure C), the harmonic content in Reference 6 is in qualitative 
agreement with experiment, considering that lower Q’s and different 


E's were observed compared to Kuntsal’s measurement [Ref. 6]. 
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VI. CONCLUSIONS 


The experimental apparatus, procedures, and computer programs 
are capable of providing data sufficiently precise to test the prediction 
of the theory in the absence and presence of geometrical perturbation. 

Because the strength parameters (STRM) and quality factors (Q) 
calculated in the experiment did not match the STRM and Q in Refer- 
ence 6, quantitative agreement could not be tested. However, the 
experimental results show that a geometrical perturbation alters the 
finite amplitude behavior of the cavity, and that the nature of these 
changes is in qualitative agreement with the predictions of the theory. 

The quantitative disagreements between the theory and experi- 
ment for the perturbed cavities may be because the resonance fre- 
quency of the non-family mode is not tracking the family modes. 

In order to avoid disturbing the position of the piston in the port 
of the cavity, the microphone has to be shifted very carefully from one 


port to another between the “pauses” as the program goes on. 
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APPENDIX A 


CURVES 


Average Values for Qs and Es 


Mode: O10 Q, = 194.72 E, = 0.0 

STRM: .389 Q, = 278.26 Eos. = 96/038 10> 
Q, = 3595.18 ea 16:36 x 10* 
Ore 205.12 ae 2.01 a0” 





Frequency Parameter 


Figure A 


Unperturbed Cavity 
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Average Values for Qs and Es 


Mode: 010 Oe = a2) E, = 0.0 

STRM: .355 Q, = 287.86 FE, = 20.39 x 10° 
Q, = 340.61 E, = 005/70 10] 
@ oo = 202.99 Eo = 23.75 x 10% 





Frequency Parameter 


Figure B 


Linear Perturbed Cavity 
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Average Values for Qs and Es 


Mode: 010 Q, = 187.52 E, = 0.0 

STRM: .370 Q, = 288.13 E, = 49.02 x 10° 
Q, = 350.42 ee 1-20.25 % 107 
Qioo = 291.40 Eo = 701.34 x 107 
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Frequency Parameter 


Figure C 


Wedge Perturbed Cavity 
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APPENDIX B 
TABLES 


TABLE 1 


UNPERTURBED CAVITY 
(SEE FIGURE A) 


Pre-—Run 


969.801 962.875 964.351 0.0 


1131.840 | 1127.791 1129.812 98.16x10-9 
1698.271 | 1693.524 1695.901 16.82x10-4 
1131.872 | 1127.590 | 1129.731 91.17x10-5 





Finite-—Amplitude Run 
V; = -11 dB Mode: 010 STRM = .389 
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of | 


es 


Post-Run 


966.381 963.270 964.82 0.0 


1132.740 | 1128.676 1130.71 94.00x10-9 
1699.770 | 1694.974 | 1697.37 16.90x10-4 
1132.800 | 1128.487 | 1130.64 74.16x10-95 
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lege)! 7 


LINEAR PERTURBATION CONFIGURATION 
(Sie PIGURE B) 


Pre-—Run 
























969.465 966.545 965.653 
11Sser8Se yu Lizeeson 132306 
1703.014 | 1698.008 1700.594 
L131s650 | 1 ladies 1134.651 


ON: 
17.672 107+ 
21 ADsaOe4 
263-021) Zoe lOee 







Finite-Amplitude Run 


Vie= a2 ds Mode: O10 SPRM = "3e5 
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970.410 967.456 966.861 





0.0 
1135.067 | 1131.116 1133.984 Da hOxanOn 
1705.321 | 1700.340 1701.110 39.98x10-4 
1136.740 | 1132.416 1135.756 17.94x10-4 


43 


TABLE 3 


WEDGE PERTURBATION CONFIGURATION 
(SEE FIGURE C) 


Pre-Run 


066.912 063.864 965.381 0.0 


1127 LSI seZen Wh Pao) Fae tO, —48.92x10-4 
1695.310 | 1690.471 1692.89] -19.17x10-4 
1Z6.76 1 al 22 ee 1124.829 -52.47x10-4 





Finite-Amplitude Run 
V; = -11 dB Mode: O10 STRM = 370 
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966.950 563.963 965.464 0.0 
Pierro |) 1123.41 1 1125.371 —49,12x10-4 
1695.544 | 1690.731 1693.144 -27.33x10-4 
1126.981 | 1123.126 1125.051 -50.21x10-4 
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TABLE 4 
E'S DIFFERENCES BETWEEN RUNS 





| contiguration | “Configuration | “Configuration 
Configuration Configuration Configuration 
vor | store | sasxtot | o2ei0-t 
pars | coseio-t | ssoxio-t | s.t6x10-4 
[Afeneta | SIO 







.53x10-4 
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APPENDIX C 
COMPUTER PROGRAMS 


This program drives the HP 3325A Synthesizer/Function genera- 
tor between 562 Hz and 1710 Hz, finds the resonance frequency, and 


calculates the quality factor (Q). 


10 | SUBROUTINE Q 

14 SENmME 0.87200 

TS FOR INC=1 TO 30 

20 K=—1 

30 FOR FREQ=562 TO 1710 STEP 2 

50 IF FREQ>550 AND FREQ<570 THEN GOTO 210 
70 IF FREQ>570 THEN GOTO 90 

90 POR FReE@=1116 TO 1140 STEP 2 

110 IF FREQ>1115 AND FREQ<1140 THEN GOTO 210 
130 IF FRQ>1140 THEN GOTO 1791 

140 K=K+1 

iO P@rrrmeo-16c0 101710 STEP 2 

170 IF FREQ>1670 AND FREQ<1710 THEN GOTO 210 
190 IF FREQ>1710 THEN GOTO 1800 

210. Fl=-RE@Q 

230) F2=Fle2 

250 F3=F1+4 

“Om REMOTE 717,723,702 

2eO) OUTEUT 717 FR Flo HZL 

310 WAIT 5000 

330 ENTER 723 :AMP1 

soo OUTPUT 717 7 FR':F2"HZ 

370 WAIT 5000 

S00 ENE R 723 > AMF2 

Bio vOUmmewT 717: FR FS HZ 

430 WAIT 5000 

450 ENTER 723; AMP3 

470 DISP F1,F2,F3,AMP1,AMP2,AMP3,K 

510 IF AMP2<AMP1 OR AMP2<AMP3 THEN GOTO 590 
570 IF AMP2>AMP1 OR AMP2>AMP3 THEN GOSUB BINOM 
SOO ee Ne eRe 

610 !SUBROUTINE BINOM 
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630 
650 
670 
S90 
710 
720 
730 
750 
770 
(328 
800 
810 
830 
850 
860 
870 
890 
ow 
930 
250 
ony 
g98 
1010 
1030 
1050 
1060 
1065 
1070 
1090 
1100 
1110 
1130 
Lee 
1170 
al SNE. 
a) 
1230 
1250 
1270 
1290 
ESO 
1330 
1350 
ei) 
1390 
1410 
1430 


BINOM: 

P=F2 

S=2 

S=S/2 
OUTPUT 717. Reece 
WAIT 5000 

ENTER 723 ;: AP 

IF P>1710 THEN GOTO 1800 
PU=P+S 
OUTRUN 717 “FR FUCA 
WAIT 5000 
ENTER 723 ; AU 

PL=P-S 
OUlPUT 17 “FR -PibanZ 
WAIT 5000 

ENTER 723; AL 

DISP P,PU,PL,AP,AU,AL 

IF AU>AP THEN P=PU 

IF AL>AP THEN P=PL 

IF S>.008 THEN GOTO 690 
GOSUB HPWR 

RETURN 

| SUBROUTINE HALFPOWER 
HPWR: 

DISP P,AP 

DISP TIME 

PRINT TIME 

PRINT "MAX" 

PRINT P,AP 

IF INC>1 AND INC<30 THEN GOTO 1771 ELSE GOTO 1110 
PRINT "REAL HALFPOWER" 
APH1=AP/24.5 

PRINT APH1 

FOR F=1 TO 2 
PH=P+((-1)F) 

OUTPUT 717 7FR':PH’HZ 
WAIT 5000 

ENTER 723; AH 

DISP PH,AH 

IF AH<APH1 THEN GOTO 1350 ELSE GOTO 1310 
PH=PH+((—1)4F) 

GOTO 1210 

S=.5 

PHI=PH-—(S*((-1)“F)) 
FOR l=1 TO 10 

OUTRUN Vi = EA ea Z 
WAIT 5000 
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1450 
1470 
1490 
Woy’ 
1530 
SSO 
1570 
gO 
1610 
1630 
1650 
1670 
1960 
ae, 
1730 
1720 
17 0 
7 
I 
773 


Tes 
Wes 
ao 2 
13s 
1794 
1 STA 
1800 
1810 


ENemer2o - An 

S=$/2 

IF AHI>APH1 THEN GOTO 1570 

PHI=PHI-(S*((-1)4F)) 

NEXT | 

GOTO 1610 

PHI=PHI+(S*((—1)F)) 

NEXT | 

PRINT "HALFPOWER" 

DISP PHI,AHI 

PRINT PHI,AHI 

IF F=1 THEN PHIR=PHI 

NEXT F 

Q=P/(PHI-PHIR) 

PRINT "Q VALUE" 

PRINT Q 

DISP Q 

K=K+1 

IF K=0 THEN GOTO 90 

Tapeh FIRS Ease omiFl THE MICROPHONE TO THE PORT CTO 
MEASURE THE FREQUENCY OF THE DEGENERATE MODE. 

PARE seCOND PAUSE SHIFT THE MICROPHONE TO THE PORT B 

IF K=1 THEN PAUSE 

IF K=2 THEN PAUSE 

IF K=1 THEN GOTO 90 

ha timenN GOTO 140 

RETURN 

NEXT INC 

END 
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This program drives the HP 3562A Dynamic SignalAnalyzer by HP 
3325A Synthesizer/Function generator and finds and computes the 
harmonic content of the microphone output at each driving frequency 


between 563 Hz and 570 Hz. 


10 REMOTE 720,717 

iS SETTIME 0,87033 

20 OUNRUT 720 4st 

30) > OURREI 20 20 

40 OUTPUT 720 ;"FRS3KH2Z" 

50°)~— COURT 720 Roe. 

60 P=800 

70 OUTPUT 717 ;"AMP";P;"MV" 

80 OUTPUT 720 ;"LNRS" 

100 @OUiRUT 72a Seer 

110. CURRED 7/208 Porc 

120° *@QUMRe7 20% Ney cee 

125 !DRIVING THE DYNAMIC SIGNAL ANALYZER AND READING THE 
HARMONIC CONTENT OF THE DRIVING FREQUENCY. 

130) Ge@-5o3 

14) VOUIE Oia? ER FREQ hz 

so SORTPWim2o 7 Ssksssnze 

160 WAIT 3000 

170 SOW Ui 20> Siri 

180 WAIT 8000 

190 OUTPUT 720 ;"MKPK" 

200 OUTPUT 720 ;"RDMk" 

210 JENNER 720. x1, VS 

220 eels 

230 BDISE Ar wile 

234 PRINT TIME 

235 PRINT X1,DIL 

240 FOR ERE@=563 TO S70isieres 

250 OUTPUT 720 ;"FRS3KH2" 

251 MOUMEU e720. Sk ERne@enz 

295 OUNEOt 7 (PR EReE@s ha 

260 FOR I|=1 TO4 

270 OUTPUT 717 ""AMP";P;"MV" 

2505 OU 20. oan 

290 WAIT 8000 

SOU) OUGEUIE 720 “MARK 


20 


310 
320 
324 
S20 
326 
Soe 
333 
335 
350 
351 
355 
356 
360 
369 
370 
380 
SN eke) 
390 
400 
410 
420 
430 
440 
450 
460 
470 
475 
480 
490 
910 
660 


OUTPUT 720 :"RDMK" 

ENTER 720 ;X,YB 

DISP X,YB 

PRINT X,YB 

PRINT TIME 

IF l=1 THEN GOTO 333 ELSE GOTO 350 
[fiKn=©>566 IhEN GOnOsce5 ELSE GOL@ 350 
GOSUB BINOM 

OuwnRPUT 720: SE’ X+50,7 AZ 

CUiirel 720; FAS. 7KAZ' 

NEXT | 

NEXT FREQ 

| SUBROUTINE BINOM 

| THIS SUBROUTINE KEEPS THE FUNDAMENTAL CONSTANT. 
BINOM: 

RAT=104(DIL/10)/104(YB/10) 

DISP RAT,DIL,YB 

IF RAT>.999 AND RAT<1.001 THEN RETURN 
P=P*RAT*%.2 

OUTPUT 717 :"AMP":P;"MV" 
SUiaeiereO “Sie 

WAIT 8000 

OUTPUT 720 ;*"MKPkK" 

OUTPUT 720 ;"*RDMK" 

ENTER 720 :X,YB 

DISP X,YB 

PRINT TIME 

PRINT X,YB 

IF FREQ>570 GOTO 660 

GOTO 380 

END 
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